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ABSTRACT 

Extended steep-spectrum radio emission in a galaxy cluster is usually associated with a re¬ 
cent merger. However, given the complex scenario of galaxy cluster mergers, many of the 
discovered sources hardly fit into the strict boundaries of a precise taxonomy. This is espe¬ 
cially true for radio phoenixes that do not have very well defined observational criteria. Radio 
phoenixes are aged radio galaxy lobes whose emission is reactivated by compression or other 
mechanisms. Here, we present the detection of a radio phoenix close to the moment of its 
formation. The source is located in Abell 1033, a peculiar galaxy cluster which underwent 
a recent merger. To support our claim, we present unpublished Westerbork Synthesis Radio 
Telescope and Chandra observations together with archival data from the Very Large Array 
and the Sloan Digital Sky Survey. We discover the presence of two sub-clusters displaced 
along the N-S direction. The two sub-clusters probably underwent a recent merger which is 
the cause of a moderately perturbed X-ray brightness distribution. A steep-spectrum extended 
radio source very close to an AGN is proposed to be a newly bom radio phoenix: the AGN 
lobes have been displaced/compressed by shocks formed during the merger event. This sce¬ 
nario explains the source location, morphology, spectral index, and brightness. Finally, we 
show evidence of a density discontinuity close to the radio phoenix and discuss the conse¬ 
quences of its presence. 

Key words: galaxies: clusters: individual: Abell 1033 - large-scale structure of Universe - 
radio continuum: general - X-rays: galaxies: clusters 


1 INTRODUCTION 

Galaxy clusters form via a hierarchical sequence of major and mi¬ 
nor mergers. A merger event dissipates up to 10^^ - 10^"^ ergs into 
heating the gas but also through large-scale intra-cluster medium 
(ICM) motions such as shocks and turbulence. These processes 
are related to the formation of non-thermal components in the 
ICM: cosmic rays (CR) and magnetic fields (e.g. |Brunetti & Joii^ 
|2014| ). Relativistic particles diffusing into magnetic fields emit syn¬ 
chrotron radiation. Therefore, non-thermal components in galaxy 
clusters are most evident once their radio emission is traced. Cur¬ 
rent observations strongly support this idea of a tight connection 
between cluster-scale synchrotron emission and mergers, finding 
indeed these radio sources always hosted in dynamically disturbed 
systems ( [Cassano et al.|2010||2013| l. 

Radio relics and radio halos are two types of radio sources 
associated with non-thermal cluster components. Radio relics (or 
gischt) are located in cluster outskirts. They are long (1-2 Mpc) 
and thin (~ 100 kpc) sources. They are polarized and are likely 
signatures of electrons accelerated by large-scale accretion shocks 
through, for example, diffusive shock acceleration 


[Blandford & Eichler] 1987 [pones & Ellison|1991| l. Radio halos are 
diffuse, unpolarized, centrally located sources. They may be caused 
by large turbulent regions formed during cluster collisions and can 
reaccelerate populations of long-lived synchrotron-dark CR elec¬ 
trons to radio bright energies ( [Donnert et al.||20T^ . These pop¬ 
ulations are generated by AGN activity, shocks or supernovae in 
galaxies ( [Brunetti & Jones|2014| l. 


A third, less studied, type of emission in merging clusters are 
the so-called radio phoenix es (|Slee et al.|2001| Ensslin & Bruggen) 


|2002[|Kempner et al.|2004[|van Weeren et al.|2009| l. When shocks 
from a cluster merger pass through an old AGN lobe, they compress 
the fossil radio plasma. Because the density in the fossil plasma 
is much lower than the surrounding ICM, these shocks may be¬ 
come subsonic compression waves inside the plasma, as shown by 
[Ensslin & Gopal-Krishna| ( [2001| ). This compression is able to re¬ 
energizing the electrons in the fossil plasma so they shine again at 
radio wavelengths. 


Here, we present the detection of a radio phoenix located in 
Abell 1033, a peculiar galaxy cluster which underwent a recent 
merger. In Sec. we describe the radio, optical and X-ray obser- 
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vations of the galaxy cluster and related data reduction. In Sec. 
we present the results of these observations and we make a com¬ 
prehensive description of Abell 1033. Discussion and conclusions 
are in Sec. and Sec. respectively. Throughout this paper we 
assume a A cold dark matter (ACDM) cosmology with Hq - 13 
km s"^ Mpc"^ Clm - 0.27 and Da = 0.73. At the redshift of Abell 
1033’s brightest cluster galaxy (z = 0.1259 + 0.0006; [Allen et al.| 
|1992j ) 1" corresponds to 2.171 kpc. All images are in the J2000 
coordinate system. 


2 OBSERVATIONS 

The study of Abell 1033 has been carried our using new Wester- 
bork Synthesis Radio Telescope (WSRT, 1.4 GHz) observations 
and archival Chandra, Very Large Array (VLA, 1.4 GHz, B- 
configuration), and Sloan Digital Sky Survey (SDSS) data. 


2.1 WSRT 

Abell 1033 was observed with the WSRlQon June 6, 2014, for 
12 h with the default 21 cm setup. Only 9 antennas participated in 
the observation because of ongoing upgrades. A total bandwidth 
of 160 MHz was recorded, spread over eight spectral windows of 
20 MHz in bandwidth, each having 64 channels. All four linear 
polarization products were recorded. The calibrators 3C147 and 
3C286 were observed at the start and end of the main observing 
run, respectively. 

We calibrated the data with CASiA[^ The first step in the 
data reduction consisted of the removal of time ranges affected 
by shadowing. We then flagged radio frequency interference us¬ 
ing the AOFlagger ( [Offringa, van de Gronde & Roerdink||2012j t. 
Before flagging, data were corrected for the bandpass response 
using 3C147. This prevents flagging of data due to the bandpass 
rolloff at the edges of the spectral windows. After flagging the data 
with AOFlagger, we again determined the bandpass response using 
3C147. We then proceeded with standard gain calibration using the 
jPerley & Butler j ( |2013j ) flux scale. The channel dependent polar¬ 
ization leakage terms were calibrated using 3C147, assuming the 
source is unpolarized. The channel dependent polarization angles 
were set using 3C286. As a last step all calibration solutions were 
transferred to the target field. 

For the target field, we carried out several cycles of phase-only 
self-calibration on 3 min timescales. This was followed by a final 
round of amplitude and phase self-calibration on 20 min timescales 
(pre-applying the previous phase-only self-calibration on shorter 
timescales). The data was imaged taking the spectral index into ac¬ 
count during the deconvolution (i.e., nterms=2 |Rau & Cornwell j 
[201 Ij ). The final image (see Fig. uses uniform weighting and 
it is corrected for the primary beam attenuation. The resolution is 
19"xl 1" and the rms noise is cr = 44 /rJy beam“^. A low-resolution 
image has been obtained after wu-subtracting the high-resolution 
clean components and tapering the data. This image has a resolu¬ 
tion of 41"x34", an rms noise of 100 /rJy beam“^and is shown in 
Fig.0 


2.2 VLA 

Searching the VLiA[^ archive we found an observation at 1425 MHz 
taken on 18th April 1993 with the “old” VLA in B-configuration 
(project code AA150). The target was observed for a total of 30 
minutes using NVSS J112553-^261020 as phase calibrator and 
3C 286 as flux calibrator. Fluxes have been scaled to match the lPer-l 
[ley & Butlerj p013j ) flux scale. Data were taken in two spectral 
windows centred on 1.39 GHz and 1.46 GHz. The bandwidth was 
of 50 MHz in both cases. After a standard calibration made with 
the CAStAP package, phases were self-calibrated and the visibili¬ 
ties were imaged using the CLEAN task with a Briggs weighting 
of -1-1 to enhance the extended emission. The final image (Fig.jlbj) 
has a resolution of 5" and an rms noise of 47 //Jy beam"^ 


2.3 Optical Imaging and Spectroscopy 

We use photometry and spectroscopic/photomet ric redshifts from 
the 10* release of the Sloan Digital Sky Survey ( |Ahn et al.|2014| 
to investigate the dynamics of Abell 1033. We limit our optical in¬ 
vestigation to a 20' X 20' (2.7 x 2.7 Mpc^) field-of-view centred 
on RA = 10^31”^42® and Dec = -h35°03'36". The brightest cluster 
galaxy (BCG) does not have an SDSS spectroscopic redshift, how¬ 
ever it was observed by [Allen et al.| ( [T992] ) which found a redshift 
of 0.1259+0.0006. This is the only galaxy within the 20'x20' field- 
of-view with a redshift in the NASA/IPAC Extragalactic Database 
(NED^that does not have an SDSS spectroscopic redshift. We in¬ 
clude this redshift with the SDSS spectroscopic redshift catalogue 
for our redshift analysis. In total there are 51 spectroscopic red- 
shifts in the field-of-view, 33 of which are within +3000 km s“^ of 
the cluster mean spectroscopic redshift (z = 0.1230+0.0005). How¬ 
ever, unless specified, in what follows the redshift of the cluster is 
assumed to be the one of the BCG (z = 0.1259 + 0.0006). We made 
this choice since there appears to be a bimodal redshift distribution 
(see Sec. |3.1.2j ) and peak of the X-ray emission and radio phoenix 
are both close to the BCG in projection. 


2.4 Chandra 

Abell 1033 was observed with Chandra for 30 ks on 2013 Eebru- 
ary 19, and for 34 ks on 2013 Eebruary 21 (PI: E. Gastaldello, 
ObsIDs 15084 and 15614). The two datasets were retrieved from 
the Chandra Data Archive, and analysed using ciao v4.5.9 and the 
Chandra Calibration Database (CALDB) v4.5.9 released on 2013 
November 19. Both observations were taken in VEAINT mode, 
with the ACIS-I and ACTS-S3 CCDs on. The level 1 event files 
were reprocessed and screened for potential bad events associated 
with cosmic rays (CR). The new level 2 event files were afterwards 
screened for background flares using the lc_clean script. Point 
sources were identified and removed from the ACIS-S3 CCD, and a 
good time intervals file was created from the source-free ACIS-S3 
event file by calculating the mean 2.5-7 keV event rate from time 
bins within 3cr of the global ACIS-S3 mean in the same energy 
band, and by then excluding all the time bins in which the event 
rate deviates from the calculated mean by a factor larger than 1.2 
( [Markevitch|2002| ) . The good time intervals defined by this screen¬ 
ing were used to remove periods of strong background flares from 


^ WSRT: http: //WWW. astron.nl ^ VLA: http://www.vla.nrao.edu 

^ CASA version 4.2.1: http://casa.nrao.edu ^ NED: http ://ned. ipac. cal tech, edu 
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Figure 1. Left: WSRT image at 1.4 GHz, contours are at [-1,1,2,4,8,16,32,64] x 3cr, with cr = 44 //Jybeam“^(beam = 19"xll"). Right: VLA image at 
1.4 GHz, contours are at [-1,1,2,4, 8,16,32,64] x 3cr, with cr = 47 juJy beam“^(beam = 5"x5"). In both images red contours trace the X-ray emission as 
described in Fig.[^ 


all the CCDs that were on during the observations. After the flare¬ 
filtering, the two observations had clean exposure times of 22.8 and 
23.6 ks. 

Count images were created in the energy band 0.5-7 keV for 
each of the two Abell 1033 datasets, and binned by a factor of 2 
(1 pixel ^ 1"). We also generated corresponding point spread func¬ 
tion (PSF) maps with an encircled counts fraction of 90%, using an 
APEC spectral model with a temperature T = 6.5 keV, a metallic- 
ity Z = 0.25 solar, and a redshift z = 0.1259; as discussed in Sec¬ 
tion |3.3.1| these parameters provide a relatively accurate descrip¬ 
tion of the global ICM parameters. Considering the same spectral 
model and energy band, exposure maps were generated for each 
observation. Using the exposure maps and the PSF maps of the 
two datasets, we created a common exposure map-weighted PSF 
map. We also created summed exposure and count images, which, 
together with the exposure map-weighted PSF map, were used to 
detect point sources in the Abell 1033 field of view with the ciao 
task wavdetect. Point sources were detected using wavelet scales 
of 1,2,4,8,16 pixels and excluded in elliptical regions scaled by 
a factor of 5cr. All point sources were confirmed visually and ex¬ 
cluded from the data, together with two additional point sources 
missed by wavdetect. West of Abell 1033, near the edge of the 
Chandra FoV, there are three galaxy cluster candidates listed in the 
NED database: GMBCG J157.77147+35.12026 (z=0.3530, |Hao 
et al.|2010J, WHL J103103.4+350623 (z = 0.3312|Wen, Han & Liu 
andNSCS J103107+350552 (z = 0.1600 |Lopes etal.|2004t . 
The locations of these cluster candidates roughly coincide, and they 
are associated with an X-ray bright AGN at RA = 10^31"^4.8® and 
DEC = 35°06'50.6". While X-ray emission from another galaxy 
cluster does not appear to contaminate the Abell 1033 FoV, to be 
conservative we excluded from the analysis a circular region with 
a radius of 3 arcmin around the bright AGN. In Fig. we show 


the exposure-corrected, vignetting-corrected Chandra map, with 
the excluded regions overlaid. 

Background event files associated with the two observations 
were created from Group “E” stowed background event files (ef¬ 
fective since 2005-09-01X00:00:00) that were reprojected to match 
the observations, and normalized such that their 10-12 keV count 
rates match the 10 -12 keV count rates recorded in the correspond¬ 
ing Abell 1033 datasets. 

The sky foreground and background emission was modelled in 
a region at the edge of the ACIS-I FoV, assuming that it is described 
by emission from the Local Hot Bubble (LHB), the Galactic Halo 
(GH), and the Cosmic X-ray Background (CXB). The foreground 
components (i.e., the LHB and the GH) were modelled with APEC 
components, while the CXB was modelled with a power-law. The 
sky background spectra extracted from ObsIDs 15084 and 15614 
were fitted in parallel, leaving the temperature and normalizations 
of all components free in the fit, but linked between the two spectra. 
The metallicities of the LHB and GH were fixed to the solar val¬ 
ues (assuming [Anders & Grevesse|1989| abundances), the redshifts 
were fixed to zero, and the power-law index of the CXB was fixed 
to 1.41 ( |De Luca & Molendi|2003j l. The hydrogen column density 
was set to the Galactic value of L64x 10^*^ atoms cm"^, as measured 
in the direction of the cluster from the Leiden/Argentine/Bonn Sur¬ 
vey of Galactic Hi ( jKalberla et al.||2005] ). To constrain the low 
temperature LHB component, we fitted simultaneously with the 
sky background spectra a ROSAT All-Sky Survey (RASS) back¬ 
ground spectrunj^ extracted from a circular region with a radius of 
1 degree around the position of the cluster centre. The CXB nor¬ 
malization of the RASS spectrum was fixed to 8.85 x 10"^ pho- 


^ http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/xraybg/ 
xraybg.pl 
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^ temperature, in units of keV; 

^ spectral normalizations, in units of pho¬ 
tons keV“^ cm“^ s“^ arcmin“^ measured at 1 keV for 
the CXB, and in units of cm“^ arcmin“^ for the LHB, 

GH, HF, and ICM; 

^ fixed parameter. 

Table 1. Best-fitting sky background parameters calculated from the two 
Chandra ObsIDs using different spectral components - first with LHB, GH, 
and CXB components, and second with an additional HF component. 


tons keV ^ 

[ 200 ^ . 


arcmin measured at 1 keV (Moretti et al. 


The spectral fits to the sky background are improved by the ad¬ 
dition of an extra APEC component with a temperature of ~ 1 keV, 
similar to the Hot Foreground (HF; |Yoshino et aLl|2009|l compo- 
nents detected towards, e.g., the Perseus Cluster {[Simionescu et aL| 
|20TT) and CIZA J2242.8+5301 jOgrean et al.|2013r However, un- 
like these other clusters, Abell 1033 is seen at a higher Galactic lat¬ 
itude. We performed an F-test to evaluate the significance of the HF 
component in the direction of Abell 1033, and found a F-test prob¬ 
ability p = 1.8% (meaning a better fit with the HF component). We 
also tried fitting the sky background spectra assuming that the pre¬ 
sumable HF emission is instead residual emission from the cluster; 
however, while this fit had an F-test probability p = 0.69% when 
compared to the fit without a HF/ICM component, we were able 
to set only an upper limit on the CXB normalization. Abell 1033 
has a virial radius /? 2 oo = 1-03 Mpc ( |Rines et al.|2013| l. The back¬ 
ground region is extracted from radii larger than ~ 0.96 Mpc from 
the cluster centre. ICM emission at ~ R 200 is not typically detected 
in ~ 45 ks Chandra observations, and therefore we consider it more 
likely that foreground emission is responsible for the hot thermal 
component detected in the sky background spectra. 

The fits were done in the energy band 0.5 - 7 keV, using Xspec 
vl2.8.2. The two sets of best-fitting parameters - without and with 
a HF component - are summarized in Table 


Figure 2. Left: Chandra surface brightness map, in the energy band 
0.5 - 7 keV. The image was binned by 2, exposure-corrected, vignetting- 
corrected, and smoothed with a two-dimensional Gaussian kernel of cr = 
3 pixels (1 pixel ^ 1 arcsec). Centre: the image shows point sources de¬ 
tected by wavdetect and confirmed visually, as well as two additional point 
sources missed by wavdetect and a circular region of radius 3 arcmin lo¬ 
cated west of the cluster, centred approximately on the position of three 
galaxy cluster candidates listed in the NED database. Right: image with the 
point sources subtracted zoomed in on the main cluster emission. 


3 Abell 1033 


The galaxy cluster Abell 1033 is at z ~ 0.12 59 and has an SZ esti- 
mated mass of M 500 = (3.4 + 0.4) x 10^"^ Mpl Planck Collaboration 


|2013| ). In this section we will present the analysis of the cluster’s 
galaxy distribution in projection and redshift. We will also study 
the radio environment and the X-ray properties of the cluster. 
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3.1 Optical properties 

3.1.1 Projected Galaxy Distribution 


We use SDSS photometry and photometric redshifts to select likely 
cluster members in order to study the projected galaxy distribution 
of the cluster. We use the SDSS type flag to disambiguate stars and 
galaxies in the held. We select all galaxies with a photometric red- 
shift within +0.103 (1 + Zciuster), which is the root-mean-square of 
the photometric redshift uncertainty based o n the SDSS faint spec¬ 
troscopic calibration sampl^ As noted in we And the mean 
cluster redshift Zciuster = 0.123. Since the photometric redshifts be¬ 
come unreliable for faint galaxies, due to photometric uncertainty 
and lack of spectroscopic training galaxies, we apply a magnitude 
cut of r < 21 to our cluster sample. Examining the w - r color- 
magnitude distribution of these galaxies we And that a number are 
redder than the two BCG’s (w - r ~ 3), thus we additionally remove 
any galaxy with u - r > 3.5. We deflne the 350 galaxies that pass 
these selection cuts as cluster members. 

We plot the distribution of cluster members, after smoothing 
with a Gaussian kernel with cr = 41", in Fig.|^ We And that the 
distribution of cluster galaxies is largely bimodal with evidence for 
a north and south subcluster. Both north and south BCGs (BCG 


north: RA 10^31”^39.8^ - DEC +35°05'34.4", BCG south: RA 
10^31"^44.3® - DEC +35°02'29.0") are within 50 kpc projected dis¬ 
tance form the relative subcluster number density peaks (subclus¬ 
ter north: RA 10^31”^40.8^ - DEC +35°05'27", subcluster south: 
RA 10^31^42.7^ - DEC +35°02'41"). The X-ray luminosity peak 
(RA 10h31”^44.8^ - DEC +35°02'18") is within 100 kpc projected 
distance from the southern subcluster number density peak, but it 
is very close to the southern BCG. Despite the similar extent and 
galaxy densities of the north and south subclusters, there is no ap¬ 
parent X-ray peak associated with the northern one. 

To examine the signiflcance of the north and south subclus¬ 
ters we created 1000 bootstrap realizations of the galaxy number 
density map. Using the variance of these bootstrap realizations and 
assuming Gaussian statistics we And that the southern subcluster 
has a peak signal-to-noise-ratio (SNR) of 5.6, the northern subclus¬ 
ter has a peak SNR of 4.9, the 150 galaxies Mpc“^ level has an 
average SNR ^4.0, and the trough region between the two peaks 
has an average SNR ^ 4.5. We also And the same structures in a 
luminosity weighted density map. Thus we believe Abell 1033 is 
composed of two subclusters of comparable richness. 


3.1.2 Redshift Analysis 

We examined the redshift properties of the 33 spectroscopic clus¬ 
ter galaxies. We show the redshift distribution of these galax¬ 
ies in Eig. 1^ Eor these galaxies we estimate a redshift of z = 
0.1230 + 0.0005 and velocity dispersion of 800 + 80kms"^ us¬ 
ing the biweight-statistic and bias-corrected 68% confldence lim¬ 
its ( [Beers, Flynn & Gebhardt|1990] l applied to 100,000 bootstrap 
samples of the cluster’s spectroscopic redshifts. Given the bimodal 
galaxy number density (Fig. and disturbed redshift distribution 
(Fig. it is likely that the velocity dispersion is biased. The red¬ 
shift distribution. Fig. shows signs of being bimodal, however 
we note that it is unlikely that the observed distribution is an accu¬ 
rate representation of the full cluster redshift distribution given the 
complex nature of the cluster and limited number of spectroscopic 
redshifts. 



10.0 10:32:00.0 50.0 40.0 30.0 20.0 31:10.0 

Right ascension 


Figure 3. Abell 1033 galaxy number density map (blue scale), based on the 
cluster member selection cuts described in |3.1.1| The blue contours be¬ 
gin at 100 galaxies Mpc“^ and increase in increments of 25 galaxies Mpc“^. 
The distribution is largely bimodal showing evidence for a north and south 
subcluster. The X-ray gas (white contours; from Fig.[^ is peaked near the 
southern subcluster, while there is no apparent gas peak near the northern 
subcluster. The VLA radio contours from Fig.^jare shown in red. 
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Figure 4. Redshift distribution of the 33 spectroscopic galaxies within 
±3000kms“^ of the mean redshift ofz = 0.1230 + 0.0005. The redshift of 
the north and south BCG’s are denoted by the arrows. While there appears 
to be a bimodal redshift distribution it is unlikely that the observed distri¬ 
bution is an accurate representation of the full cluster redshift distribution 
given the complex nature of the cluster and limited number of spectroscopic 
redshifts. 


6 


http: //WWW .sdss3.org/dr10/algorithms/photo-z.php 
































6 F. de Gasperin et al. 


To investigate whether there is a correlation between the bi- 
modal structure in the projected galaxy number density (Fig.[^ and 
redshift distri bution (Fig.[^ we perfo rm both a Dressler-Shectman 
test (DS-test; [Dressier & Shectman||1988| ) and a simple redshift 
analysis of the spectroscopic galaxies near each subcluster. We 
perform the DS-test using { V331 nearest neighbours and find no 
evidence of subclustering. For the simple redshift analysis of the 
north and south subclusters we select the spectroscopic galaxies 
within a radius of 160" (~350kpc) of each respective subcluster 
number density peak. These circular apertures are as large as pos¬ 
sible while maintaining mutual exclusivity of the subcluster spec¬ 
troscopic members. These selections resulted in 9 galaxies being 
assigned to each subcluster. We perform the same redshift analysis 
as done for the system and find the redshifts of the north and south 
subclusters to be 0.1240 ^qqqq 9 and 0.1224^QQQQg, respectively, and 
the velocity dispersion^ to be 720^}^^ km s"^ and 690^^^^ km s"^ 
respectively. We find that the calculated redshift of each subclus¬ 
ter is located between the peaks of the redshift distribution (Fig.|^ 
suggesting that the redshift estimate of each subcluster may be bi¬ 
ased due to the mixing of two redshift populations. 

Comparing the relative line-of-sight velocitie^of the BCG’s 
and subclusters provides further evidence that there are two over¬ 
lapping redshift distributions. We find a large relative velocity dif¬ 
ference between the north and south BCG’s, t^north,BCG - i^south,BCG = 
-1776 km s"^ however we find a much smaller relative velocity 
difference, and of the opposite sign, between the north and south 
subclusters, I’north.subciuster “ ^south,subciuster — 500 + 300 km s . Simi¬ 
larly we find large relative velocities between each subcluster and 
its respective BCG, i^nortb,subcluster - i^nortb,BCG = 1300 + 200 km s"^ 
and iisoutb,subcluster “ i^soutb,BCG = “900 + 200 km s"^ Even if this is 
a post merger cluster it seems unlikely that BCG’s could receive a 
large enough dynamical kick to have such large line-of-sight veloc¬ 
ities relative to their subclusters while remaining near the subcluster 
centres in projected space. A more plausible explanation is that the 
redshift estimate of each subcluster is biased significantly due to 
mixing of two subcluster populations with large relative velocities. 
Without a more extensive redshift survey we believe that the BCG’s 
provide the least biased estimate of the two subcluster redshifts and 
relative velocities. It is worth noting that the northern BCG has a 
redshift offset from one of the peaks of the observed redshift dis¬ 
tributions (Fig.|^, but this may be due to the observed distribution 
inaccurately representing the full cluster redshift distribution. 


3.2 Radio properties 

3.2.1 Unclassified extended radio emission 

Extended radio emission is present 35" below the cluster centre 
(defined as the peak of the X-ray emission), which corresponds to 
~ 76 kpc (see Fig.[^. We named this source Source 1. The emis¬ 
sion extends 281 kpc in length and 39 kpc in width and it does not 
have any clear optical counterpart (see Fig.[^. A plume extends 
from the east side of the source towards south-east. The southern 
rim of Source 1 terminates quite sharply while the northern one 
shows an irregular border with faint extensions elongating towards 
the cluster’s centre. Using WSRT and VIA data, we investigated 


^ Again, we caution that both the redshift estimates and velocity dispersion 
estimates of the subclusters are likely to be biased and, as discussed below, 
this bias is expected to be large. 

^ We use a relativistic velocity convention (see e.g., 


Dawson 


2013 


Length 

Width 

Distance from 

cluster centre 

Flux 

Density 

Power 

(kpc) 

(kpc) 

(kpc) 

(mly) 

(WHz-‘) 

281 

39 

76 

48 + 0.5 

1.85 X 


Table 2. Details of the extended source {Source 1). Flux density and power 
are estimated at 1.425 GHz. 


Telescope 

Frequency 

(MHz) 

Flux 

(mly) 

Reference 


VLA 

1465 

45.8 + 1.3 

this work 


WSRT 

1422 

46.9 ± 7.6 

this work 


VLA 

1385 

51.2+1.5 

this work 


WSRT 

1341 

53.9 + 7.3 

this work 


WSRT 

608 

220 

Roland et al. 

1985 

Texas interferometer 

365 

380 

Roland H^al. 

1985 


Table 3. Flux densities of the extended sources {Source 1). 


the presence of polarized emission from Source 1, and found the 
radiation <1% polarized in the brightest regions. 

The flux densities for Source 1 extracted from the WSRT 
and the VLA maps are in agreement in the overlapping frequency 
range, which means that the VLA observation was able to pick 
up all the flux although the source was extended and the array 
quite sparse (B-configuration). The two observations gave a flux of 
~ 48 + 0.5 mJybeam"^ at 1.4 GHz, which corresponds to a power 
of Fi 4 oo = 1.85 X 10^"^ W Hz"^ The properties of Source 1 are sum¬ 
marized in Table while flux densities at different frequencies are 
in Table [3] 

This extended radio source was already noticed by |Roland| 
|et al.|(T985] ) while searching for steep radio sources in galaxy clus¬ 
ters. They found the spectral index of the source to be rather steep 
{a < -1.3, Fy oc y®), but the authors could not clearly classify it 
because of a missing optical counterpart. Using data from the lit¬ 
erature (see Table we extracted a global spectral index value for 
Source 1 of a - -1.62 + 0.09 (see Fig.[^. Data points were taken 
from different telescopes covering different wi;-ranges, which may 
bias some measurements. However, the spectrum shows a hint of 
curvature becoming steeper at higher-frequencies, as expected by 
standard synchrotron ageing mechanisms. 

To study the presence of a radio halo we examined the low- 
resolution WSRT map after subtracting the emission of the high- 
resolution map. The result, displayed in Fig.[^ shows an excess of 
emission throughout the cluster which is close to the noise level. 
This can be related to the presence of a radio halo, which would 
not be surprising in a system which underwent a recent merger. 
However, the emission is too faint to derive firm conclusions. 


3.2.2 Other radio sources 

The radio environment of Abell 1033 is quite complex. A head- 
tail (H-T) radio galaxy (RG), oriented north-south, is present close 
to the cluster centre. It has a clear optical counterpart 2MASX 
110314175-^3502411 {z = 0.120) which is probably part of the 
northern sub-cluster of Abell 1033. The source has been labelled 
Source 2 (see Fig.[^. 

A wide angle tail (WAT) RG is present on the east side 
of the cluster and it is labelled Source 3 (see Fig. |^. An op¬ 
tical counterpart for this object has been identified as 2MASX 
110320162-^3502530 {z - 0.121) which is also part of the north¬ 
ern sub-cluster of Abell 1033. The emission from the WAT RG is 
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Figure 5. Top-left panel: Chandra surface brightness map, in the energy band 0.5 - 4 keV. The image was binned by 4, exposure- and vignetting-corrected, 
and smoothed with a two-dimensional Gaussian kernel of cr = 2 pixels (1 pixel = 1.97 arcsec). Points sources were subtracted from the image, and each 
resulting gap was filled with values drawn from a Poisson distribution with a mean equal to that of the pixel values in a small elliptical annulus surrounding 
the gap. Because this image is highly-processed, it has not been used in the scientific interpretation. Contours are at [0.1,0.25,0.5,1,2,3, 8] x 10“^ s“^ cm“^. 
Red contours are from the VLA observation as described in Fig.[T] Other panels: Sloan digital sky survey composite image (/, r, and g filters). Blue contours 
trace the projected galaxy density distribution as in Fig.|^ White contours trace X-ray emission as shown in the top-left panel. Red contours are from the VLA 
observation (Fig, [lb) at [1,2,4,8,16,32,64] x 3cr, with cr = 47 //Jybeam“^ In the last image contours at [1.5,2.5,3.5] x 3cr are also plotted. The optical 
counterpart of the head-tail and the wide angle tail radio galaxies are detected, no optical counterpart for the extended emission (Source 1) is instead detected. 


polarized up to 15%. Fig. shows the orientation of the E-vectors 
(perpendicular to the magnetic fields) for this source. Using the 
map provided by |Oppermann et aL| ( |2012| ) we estimated the Fara¬ 
day depth in the direction of Source 5 to be ~ 9.3 rad m"^. The 


polarization angle has therefore been corrected for Galactic Fara¬ 
day rotation by subtracting an angle of 25°. 

A fourth radio source is interesting for further discussion. It is 
a faint {S 1400 = 0.5 mJy) point source (PS) with an optical counter¬ 
part (SDSS J103143.69-F350154.5, z = 0.125) located in projection 
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Figure 6. Radio spectrum for the extended source (Source 1) in Abell 1033. 
Data points are listed in Tablej^ Blue line is a linear fit. The errors of the two 
lowest frequency data-points are unknown. In these two cases, for htting 
purposes, we assumed conservative error equal to 10% of the flux density 
values. 


Surface brightness [mJy beam D 
0.2 0.3 0.4 0.5 0.7 1.0 



50.00s 40.00s 

RA 02000) 


10h31m30.00s 


Figure 7. WSRT low-resolution map of the Abell 1033 galaxy cluster. 
The high-resolution image (shown in red contours. Fig. im has been uv- 
subtracted and data were tapered to the resolution of 41"x34" to en¬ 
hance the residual extended emission. Contours are at [1,2,3,4] x 3cr with 
cr = 100 juJybeam“^ Gray contours trace the X-ray emission as in Fig.[^ 


Figure 8. WSRT 1.4 GHz map of the WAT source NVSS J103200-h350249. 
Black contours are at [1,2,4,8,16,32,64] x 3cr, with cr = 44 //Jy beam“^ 
red contours are from the VLA 1.4 GHz observation at [1,2,4,8,16,32] x 
3cr, with cr = 47 //Jybeam“^ E-Vectors are displayed with length propor¬ 
tional to the fractional polarization which ranges from 2 to 15%. A correc¬ 
tion for Galactic Faraday rotation has been applied following |Oppermann| 
|etaLH2^ . 


29 kpc North of the brightest peak of Source 1. Given its location 
and redshift the source is probably part of southern sub-cluster of 
Abell 1033 and it is labelled Source 4 (see Fig. [^. A list of the 
proprieties of the sources discussed here is given in Table|^ 

3.3 X-ray properties 

3.3.1 X-ray global parameters 

To determine the global cluster parameters, we extracted cluster 
spectra from a circular region of radius 5.3' (^ 690 kpc) centred at 
RA = 10^31”^44.4^ and DEC = -f35°03'12.0". This region corre¬ 
sponds approximately to the cluster’s R 500 , and encloses the main 
part of the visible X-ray emission from the ICM. From the total 
spectra, we subtracted instrumental background spectra extracted 
from the same regions. The sky background parameters were kept 
fixed, in turn, to each of the two best-fitting models summarized 
in Table Emission from the ICM was modelled with a single¬ 
temperature APEC model with free temperature, metallicity, and 
normalization, and with a fixed redshift of 0.1259. The spectra from 
ObsIDs 15084 and 15614 were fitted simultaneously, under the as¬ 
sumption that they are described by the same model parameters. 
The fits were done in the 0.5 - 7 keV energy band. 

For the best-fitting sky background model without a hot fore¬ 
ground component, we found an average cluster temperature of 
6.34 + 0.16 keV and a metallicity of 0.26 + 0.04 solar. For the best¬ 
fitting sky background model that included a HE component, we 
found very similar ICM parameters: an average cluster temperaure 
of 6.39 + 0.16 keV, and a metallicity of 0.25 + 0.04 solar. The cluster 
has a 0.5 - 2 keV flux of 3.86 ^qq3x 10“^^ erg cm“^ s"^ The calcu¬ 
lated flux is essentially the same for both sky background mod- 


























Abell 1033: birth of a radio phoenix 


9 


Source name RA^(J2000) Dec^(J2000) ‘^ 1.425 ghz (mJy) Note 


FIRST J103143.9+350141 
FIRST J103141.7+350231 
NVSS J103200+350249 
SDSS J103143.69+350154.5 


10^31”'43^ 

10 ^ 31 “' 41 . 77 ‘ 

10h32“'01.6F 

10h31“'43.69‘ 


+35°01'42" 

+35°02'41.42" 

+35°02'52.94" 

+35°01'54.50" 


48.0 + 0.5 

46.7 + 0.2 

46.8 + 0.5 
0.5 + 0.05 


- Extended source (Source 1) 

0.120 Head-tail RG (Source2 ) 

0.121 WAT RG (Source 3) 

0.125 Point source (Source 4) 


^ coordinates are from the SDSS optical counterpart. In the case of the extended source, they are at the peak of the radio emission. 
^ spectroscopic redshift from the SDSS optical counterpart. 


Table 4. Coordinates and fluxes of the radio sources. 


Sky background 
(photons s“^ cm“^ arcmin“^) 




ai 


(arcmin) 


5.92 X 10-6 
6.60 X 10-6 
7.30 X 10-6 


0 93+^-^^ 

0 94 + 0-17 

n T Q 


1.63 + 0.12 
1.67 + 0.13 

1 79+0.14 
-0.13 


9 9 C+O.O 5 
^•^'^-O.OS 

2.25+^-^^ 
2.25' 


15 
-0.03 


^ http://wwwasdoc.web.cern.ch/wwwasdoc/minuit/nodel4.html 


FCN^ 


1 42 + 0.18 
-0.16 
1 41+0.18 
^•^^-0.16 
1 40+018 
-0.16 


-1548.07 

-1547.99 

-1547.88 


Table 5. Best-fitting broken power-law density models to the surface brightness profile extracted from the sector shown in Fig.|^ The sky background value 


was kept fixed in the fit. We ran the fit for the best-fitting sky background level of 6.60 x 10 
the boundaries of the 90% confidence interval. 


photons s 


, and for sky background levels at 



32m00.00s 54.00s 48.00s 42.00s 36.00s 10h31m30.00s 

RA (J2000) 

Figure 9. Chandra surface brightness map in the energy band 0.5 - 3 keV. 
The image was binned by 4, exposure-corrected, vignetting-corrected, and 
smoothed with a two-dimensional Gaussian kernel of cr = 2 pixels (1 pixel 
= 1.97 arcsec). Point sources were subtracted from the image, and each re¬ 
sulting gap was filled with values drawn from a Poisson distribution with a 
mean equal to that of the pixel values in a small elliptical annulus surround¬ 
ing the gap. This is shown for clarity, but because it is highly processed it 
was not used in the analysis. The overlaid sector shows the region in which 
a density discontinuity was detected by eye. The discontinuity was mod¬ 
elled with a broken power-law density model, as described in the text and 
shown in Fig.[T^ Radio contours based on the 1.4 GHz VLA radio image 
are drawn at [1,4,16,64] x 100 //Jybeam-^ 


els, and it implies that Abell 1033 has a 0.5 - 2 keV luminos ity of 
(1.48+0.01)xl0"^'^ erg s-^ Using the Msoq-^ relation found by |Pratt| 
et aL| ( |2009t we can infer a cluster mass of (3.5 + 0 . 2 ) x 10 ^"^ M©. 


3.3.2 Density discontinuities in the ICM 

The main aim of the Chandra analysis was to identify density and 
temperature discontinuities in the ICM, in order to explore their 


connection with radio sources. The only clear density discontinuity 
that we were able to detect by eye is located SE of the cluster, in 
the elliptical sector shown in Fig.[^ We extracted a surface bright¬ 
ness profile from this elliptical sector, using images and exposure 
maps in the energy band 0.5 - 3 keV. The instrumental background 
surface brightness profile extracted from the same sector and in the 
same energy band was subtracted from the total surface brightness 
profile, and binned to a minimum of 1 count/bin. We modelled the 
sky background by fitting a constant to the outer regions (13 - 18 
arcmin) of the profile. Once the sky background level was deter¬ 
mined, we fitted a broken power-law density model to the surface 
brightness profile (in the radius range 1-5 arcmin), keeping the 
sky background level fixed and assuming isothermal plasma. The 
broken power-law density model is defined as: 


n 2 {r) = C fZo — , for r <r^ 

\^d/ 


n\{r) 


no ^ —j , for r > 


( 1 ) 


where n is the electron number density, C is the density compres¬ 
sion, r is the distance from the centre, is the distance of the den¬ 
sity discontinuity from the centre, a is the power-law index, and 
subscripts 1 and 2 refer to the outer side of the density disconti¬ 
nuity and inner side of the density discontinuity, respectively. The 
density profile is integrated along the line of sight under the as¬ 
sumption that the cluster’s 3D geometry is a prolate spheroid. If 
the cluster would instead be longer along the line of sight, then the 
density jump we detect would be stronger. 

We also varied the sky background level within its 90% confi¬ 
dence limits, and repeated the fitting. All fits were done using prof- 
fit vl.2 ( [Eckert, Molendi & Paltani|2011| | and Cash statistics ( jCashj 
[T97^ . 

For the sector shown in Fig.l^ the sky background level in the 
energy band 0.5-3 keV is 6 . 60 ^q 4 ^x 10"6 photon s"^ cm"^ arcmin“^ 
(errors are quoted at the Icr level; the 90% confidence range is 
[5.92 - 7.30] X 10"6 photon s"^ cm"^ arcmin"^). The best-fitting 
broken power-law models for different sky background values are 
summarized in Table The background and ICM fits to the sur¬ 
face brightness profiles are shown in Fig.[^ The best-fitting bro¬ 
ken power-law model parameters are essentially unchanged as the 
sky background level is varied. The best-fitting broken power-law 

















10 F. de Gasperin et al. 


model derived for the best-fitting sky background level has a den¬ 
sity compression C = 1.41 which implies a very weak Mach 
number of 1.28 ^q }^ if the discontinuity is a shock front. The po¬ 
sition of the density discontinuity is shown in Fig. The bro¬ 
ken power-law fit was compared to a kink model fit (i.e. a broken 
power-law model with a fixed density compression C - 1) and to 
a simple power-law fit. Based on the likelihood-ratio test, a kink in 
the profile is confirmed at 4.2cr (confidence level 99.997%), while a 
density compression C > 1 is confirmed at 1.82cr (confidence level 
93.197%). 

We searched for other density discontinuities in the ICM us¬ 
ing narrow (~ 20- 30 degree) circular sectors that approximately 
follow the X-ray isophotes. However, no other discontinuities were 
detected with a confidence level > 90%. 


3.3.3 Spectral analysis near the density discontinuity 

To identify the nature of the density discontinuity detected ^ 2'25 
SE of the cluster centre, we extracted spectra in two partial ellipti¬ 
cal annuli on both sides of the discontinuity. The two annuli have 
widths of approximately 0.7 arcmin (inner region) and 1.1 arcmin 
(outer region). A negative temperature jump would indicate that the 
discontinuity corresponds to a shock front. Instead, if the temper¬ 
ature increases significantly on the outside of the density disconti¬ 
nuity, then the discontinuity would correspond to a cold front. 

The two spectra were binned to a minimum of 25 counts/bin. 
To fit the spectra, we subtracted the instrumental background spec¬ 
tra and fixed the sky background model to the best-fitting mod¬ 
els summarized in Table The ICM was modelled with a single¬ 
temperature APEC model with a fixed metallicity of 0.25 solar and 
a fixed redshift of 0.1259. The hydrogen column density was fixed 
to the Galactic value of 1.64 x 10^*^ atoms cm"^. The spectra from 
ObsID 15084 and 15614 were fitted simultaneously, under the as¬ 
sumption that they are described by the same physical model. The 
results of the fits are summarized in Tabled 

No temperature jump is detected across the density disconti¬ 
nuity. However, the high temperatures of the ICM and the relatively 
short exposure time of the observations make the measurements too 
uncertain to exclude a temperature jump corresponding to a Mach 
number of ~ 1.3. 


3.3.4 X-ray morphological estimators 

In order to classify Abell 1033 as a possible merging cluster we 
analysed the X-ray surface brightness inside an aperture radius of 
Rap = 500 kpc (230") using the methods described in|Cassano et al.| 

We calculated the cluster emission centroid shift and the sur¬ 
face brightness concentration parameter. The centroid shift, w, is 
computed in a series of circular apertures centred on the cluster 
X-ray peak and is defined as the standard deviation of the pro¬ 
jected separation between the peak and the centroid in unit of Rap, 
as ( [Poole et al.|2006[[Maughan et al.|2008| ): 

where A/ is the distance between the X-ray peak and the centroid 
of the iih aperture. 


[ht!] 





Figure 10. Left: surface brightness profile extracted from the sector shown 
in Fig. For clarity, the profile was binned to a uniform signal-to-noise 
ratio of 6. The blue line shows the instrumental background profile; in the 
plot, this is subtracted from the total profile. Middle: constant model fit to 
the sky background. For clarity, the profile was binned to approximately 30 
counts/bin. The fitted model is shown in blue. The bottom panel shows the 
fit residuals. Right: broken power-law density model fit to the radial range 
1-5 arcmin. For clarity, the profile was binned to a uniform signal-to-noise 
ratio of 6. The fitted model is shown in blue. The bottom panel shows the 
fit residuals. 
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Sky background T 2 N 2 ;y^/d.o.f. Ti Ni ^^/d.o.f. 

model (keV) (cm“^ arcrnin"^) (keV) (cm“^ arcmin"^) 

withHF (3.95 ± 0.15) x lO"'* 0.99 5.57+J ® b.dl+O-^^ x IQ-^ 0.66 

without HF 6.13:|;*'| (3.95 ± 0.15) x lO"'* 0.99 5.42+['^^ 6.41+™ x IQ-^ 0.67 

Table 6. Best-fitting APEC models to the spectra extracted from both sides of the SE density discontinuity. The hts were run with the sky background model 
fixed, in turn, to each of the two best-fitting models summarized in Table^ 


The concentration parameter { [Santos et al.|2008| ), c, is defined 
as the ratio of the peak over the ambient surface brightness, 5, as: 

^ 5(r<100kpc) 

5(r<500kpc) ^ 

Abell 1033 has a centroid shift w = 0.086 + 0.006 and a 
concentration parameter c = 0.200 + 0.004. These results put the 
cluster in an intermediate position being classified as disturbed for 
the centroid shift and as moderately relaxed for the concentration 
parameter ( jCassano et al.|201^ . Abell 1033 is therefore different 
from clusters undergoing major mergers. This can be caused by 
the merger happening close to the line of sight (as suggested by 
the BCGs high differential velocity) or by the mass of the merging 
clusters being very different (minor merger). A third possibility is 
that the age of the merger is comparable to the sound crossing time 
of the X-ray gas (old merger). Thus the gas might have relaxed to 
some degree. This relaxation time can be less than the time it takes 
for the two subcluster dark matter halos and galaxies to reach their 
maximum separation and come back together. 


4 DISCUSSION 

The X-ray brightness distribution of Abell 1033, together with the 
high radial velocity of the BCGs, suggests a merger axis with a 
large line-of-sight component. The compression of the X-ray con¬ 
tours (see Fig. towards south shows that the merger axis must 
also have a north-south component, which is in line with the loca¬ 
tion of the subclusters detected by optical observations. The X-ray 
emission is mostly concentrated on the south sub-clusters. An ex¬ 
planation could be that the dense south core is being ram-pressured 
stripped as it moves through the ICM, while the north cluster did 
not enter this merger with a dense core, or perhaps its core was 
badly disturbed. 

Extended radio emission {Source 7, in Fig. with a steep 
spectrum (a ^ -1.62 + 0.09) and no optical counterpart is present 
close to the region where the X-ray surface brightness gradient is 
steepest. Chandra data do not show any sign of a shock presence or 
temperature jump along what, at a first sight, could be interpreted 
as a radio relic. The extended radio source has some other pecu¬ 
liar characteristics compared to known radio relics: (i) the source 
is small (280 kpc), (ii) it has a high surface brightness, (iii) it is 
centrally-located in the cluster (unexpected for radio relics; pWzzaj 
jet al.|201^ and (iv) its radiation is not polarized down to the 1 % 
level. Given these characteristics either the source is a very unusual 
radio relic or its nature is different. 

The radio morphology is in-line with the idea of Source 1 
moving north to south. The radio brightness has a quite steep gra¬ 
dient and terminates sharply on the south edge, while the northern 
edge is disrupted and shows faint structures elongating towards the 
cluster centre. Interestingly, we discovered the presence of a den¬ 
sity discontinuity in the south-west region of the cluster beyond the 
radius of the southern edge of Source 1. The feature is located at the 



RA (J2000) 


Figure 11. Same as Fig.|^but zoomed in on the density discontinuity and 
plume region. The black dashed line shows the discontinuity location. Blue 
contours are from the low-resolution WSRT data (Fig.fl^, while red con¬ 
tours are from the high-resolution VLA data (Fig. 


edge of a radio plume, extending from the east side of the source 
(see Fig. [TT). The plume could be therefore the result of a shock- 
driven r^-acceleration of CR electrons originated from Source 1 and 
points towards an association between Source 1 and the disconti¬ 
nuity. Reacceleration of CR electrons, in contrast with direct accel¬ 
eration from the thermal pool, is thought to be required to power 
extended radio sources in galaxy clusters (e.g.jEnss lin et al.| 1998[ 

Markevitch et al.|2005||Kang & Ryu|201 l||Pinzke, Oh & Pfrommerj 

2013r 

Given the difficulties to interpret Source 7 as a radio relic, we 
present an alternative explanation for the formation of this source 
which is related to the radio emission from another cluster member 
(SDSS 1103143.69-^350154.5, Source 4 in Fig. [^, located just 
above (29 kpc in projection) the centre of the extended emission. 
The radio source, probably related to nuclear activity, is coincident 
with an early-type galaxy whose redshift (z = 0.125) and proximity 
to the southern BCG suggest it is part of the southern sub-cluster 
(see Fig. [^). The extended radio emission has a brightness and 
a morphological configuration compatible with a Fanaroff-Riley I 
radio galaxy. Therefore, it could be a striking example of CRs pro¬ 
duced by an old AGN outburst, displaced, and likely compressed 
by the galaxy cluster merger event through shock waves (i.e. a ra¬ 
dio phoenix). While the thermal environment gets shocked, the ra¬ 
dio plasma itself is only compressed adiabatically due to the much 
higher internal sound speed. Plasma compressed to become a radio 
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Figure 12. Possible configuration of Abell 1033 view from the side. The 
phoenix extensions is perpendicular to the line-of-sight. The observer view¬ 
point is on the right. 


phoenix close to the cluster centre is expected to produce bright 
emission with a very steep spectrum, as observed in this case ( |En-| 
[sslin & Gopal-Krishna|2001| l. This scenario would explain the non¬ 
detection of a shock along its extension, the source location, mor¬ 
phology, spectral index, and brightness. High degrees of polariza¬ 
tions are expected as the magnetic field is aligned and compressed. 
The absence of polarization may be, in this case, an observational 
limit due to the poor resolution in combination with the strong Fara¬ 
day depolarization due to the large quantity of intervening matter. 

Another explanation could be that Source 1 traces the fading 
lobes of a radio galaxy {Source 4) which became quiescent while 
moving north (similar to the case of 3C 338; |Gentile et al.|2007| l. 
However, the north and south edges of Source 1 show a movement 
of the lobes (compression south, rarefaction north) which is incom¬ 
patible with old lobes left behind by an AGN moving north. More¬ 
over, Source 1 would have had a narrow/wide angle tail morphol¬ 
ogy bended south, which is not seen in an environment where the 
other large radio galaxies show signs of interaction with the ICM. 

Another possible explanation is that Source 1 was not dis¬ 
placed by a merging shock but by the bulk ICM motions induced 
by the merger event. In the simplest scenario, when the two sub¬ 
clusters merge, their dark matter and galaxy components pass each 
other relatively undisturbed, while the ICM interacts and slows 
down creating a relative velocity between the ICM and the galaxies. 
Assuming the sub-clusters just passed each other, as suggested by 
the presence of a brightness discontinuity and by the N^S orienta¬ 
tion of the head-tail which is in the north sub-cluster (see Fig.p^), 
we would expect the ICM in the south sub-cluster to be moving 
towards the north with respect to the cluster galaxies. Alternatively 
said, the member galaxies of the southern sub-cluster are moving 
south faster than the ICM is. This would have shifted the lobes to¬ 
wards the north rather than the south of the AGN, disfavouring this 
interpretation. 

The most simple explanation that takes the above evidence 
into account is that the radio lobes have been displaced and com¬ 
pressed by a low Mach number merger shock moving outward in 
the ICM. The adiabatic compression boosts the lobes fiux as in the 
scenario A (radio cocoon close to the cluster centre) in |Ensslin~^ 
|Gopal-Krishna| ( [200T] ) . Such a revived radio lobe cannot be very old 


(< 100 Myr), thus must be close to the galaxy it originated from but 
displaced in the shock wave direction of propagation (i.e. south). 
Here we argue that Source 4 can be that source and that the merger 
has revived the radio cocoon. The detected discontinuity in the X- 
ray emission may well be (part of) the shock that originated the 
radio phoenix. Radio observations at low frequencies will provide 
a spectrum that can be used for deriving the age and the compres¬ 
sion factor (subject to a number of assumptions) confirming this 
scenario. 

Here we presented the simplest scenario which takes all the 
observable facts into account. Complications may rise if for exam¬ 
ple the ICM has relevant local velocities. Source 4 has an unusual 
peculiar velocity, projection effects mask the real shape of Source 
1, or the merger is more complex than a well-behaving dissociative 
merger (as a minor merger or a merger with a high impact parame¬ 
ter). 


5 CONCLUSIONS 

We presented Chandra, WSRT, and VLA observations together 
with the analysis of SDSS images of the galaxy cluster Abell 1033. 
The main findings of the paper are: 

• The galaxy number density and redshift analysis performed 
on SDSS images shows a clear bimodal distribution. Abell 1033 is 
composed of two subclusters of comparable richness. 

• X-ray brightness maps show a disturbed cluster according to 
global morphological parameters. We report a steepening in the sur¬ 
face brightness profile towards south. Abell 1033’s merger axis is 
almost oriented along the line of sight. This projection effect could 
be responsible for the non-detection of merging shocks. 

• The radio environment of Abell 1033 is quite heterogeneous. 
We report the presence of a head-tail radio galaxy and of a wide- 
angle tail radio galaxy which are probably part of the northern sub¬ 
cluster. The N^S orientation of the head-tail radio galaxy suggest 
a motion in line with the merger scenario. 

• An extended, steep-spectrum {a ^ 1.62 + 0.09) radio source is 
reported. The source is likely a newborn radio phoenix originated 
by a merger-induced shock-waves which displaced (and possible 
compressed) the lobes of a radio galaxy located 29 kpc (in projec¬ 
tion) above the source. 

• We discovered a density discontinuity which could be caused 
by a weak shock with a Mach number of 1.28 ^q}p probably related 
to the recent merger. Interestingly, the discontinuity seems related 
not to the radio phoenix itself, but rather to a plume extended from 
the radio phoenix (see Fig.[^. We speculate that the discontinuity- 
related shock may have re-accelerated CR electrons present in the 
old radio galaxy lobe. On larger scales, this mechanism is thought 
to be responsible for generating giant radio halos. In some cases 
giant radio halos are in fact found to be spatially connected to a 
shock on their edge (see e.g. |van Weeren et al.|20I2[|de Gasperm| 
|et al.|2014] ). 

Optical and multi-frequency radio observations of Abell 1033 
are planned and will shed light on this complex system. If the 
merger axis is slightly tilted N-S, then possible radio relics are 
expected to the north and/or to the south of the cluster. The pic¬ 
ture will be also reinforced by the detection of a radio halo at low- 
frequencies. 
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